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ABSTRACT 


In Part { (Schwiderski, 1978a) of this report, a unique hydrodynamical interpolation tech- 
nique was introduced, extensively tested, and evaluated in order to compute partial global ocean 
tides in great detail and with a high degree of accuracy. This novel method has been applied to 
construct the diurnal luni-solar declination («, ) ocean tide with a relative accuracy of better than 
5 cm anywhere in the open oceans. The resulting tidal amplitudes and phases are tabulated on a 
1° x 1° grid system in an atlas of 42° x 71° overlapping charts covering the whole oceanic globe. A 
corresponding atlas of global corange and cotidal maps is included to provide the reader with a 
quick general overview of the major tidal phenomena. The specifying hydrodynamical parameters 
of the model are listed along with quoted sources of empirical tide data, and significant tidal 
features are explained and discussed. The diurnal K, ocean tide is found to resemble qualitatively 
the semidiurnal M, and S, tides presented in Parts I],and HI of this report. However. major 
shifts of the positions of the amphidromes are apparent. 
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I. INTRODUCTION 


Part I of this report (Schwiderski, 1978a) introduced a unique combination of hydrody- 
namical and empirical methods to model detailed ocean tides with a relative component accuracy 
of better than 5 cm anywhere in the open oceans. This enormous accuracy is well above minimum 
requirements set by, for instance, the National Aeronautics and Space Administration (NASA) 
and the Department of Defense (DoD) — to map the geoid at sea by satellite altimetry to within 
10 cm. The following features of this unique hydrodynamical interpolation mode! made the 
achievement of this accuracy possible. 


a. A spherically graded 1° x 1° grid system is set up in connection with a corresponding 
1° x 1° bathymetry to assure a sufficient resolution of all important tidal phenomena. 


b. The bathymetry of the gridwise, simply connected ocean basin is hydrodynamically 
defined (Schwiderski, 1978c) by appropriate modifications of earlier realistic depth data col- 
lections. The hydrodynamical redefinition was needed in order to model the well-known strong 
distortion and retardation effects of shallow continental shelves. narrow ocean ridges or island 
chains, and other significant bottom irregularities. 


c. The Boussinesq substitution of the turbulent Reynolds stresses is applied in, the form of 
eddy dissipation with a novel physically meaningful eddy viscosity that depends linearly on the 
lateral grid-cell area and, hence, directly on the ocean depth. 


d. The linear law of bottom friction is introduced with a bottom-friction coefficient depend- 
ing linearly on the bottom grid-cell area which is independent of the ocean depth. In boundary 
cells, the otherwise constant friction coefficient is subjected to an indirect cellwise adjustment in 
order to permit a consistent hydrodynamical interpolation (see h., below) of empirical tide data 
known from tide gauge stations at continental shores, islands. or other shallow-ocean bottom 
irregularities. 


e. The effects of the terrestrial tide and the oceanic tidal load ure included ds simple second- 
order approximations in the sense of Love and Accad and Pekeris (1978). 


f. The Hansen-Zahel (Zahel, 1970 and 1977; Estes. 1977) finite differencing technique is 
modified by a new differencing scheme in time which improved decay, dispersion, and stability 
characteristics of the numerical procedure and facilitates the simple indirect adjustment of the 
bottom-friction coefficient in the hydrodynamical interpolation technique (see d. and h.). 


g. At land-ocean cell walls, the conditions of no-flow across and free-slip along the bounda- 
ries are enforced. The no-flow condition is subsequently relaxed by allowing controlled periodic 
inflows and outflows over the mathematically assumed boundaries. This allowance redefines 
indirectly more realistic shorelines in order to further improve the consistency of the hydro- 
dynamical interpolation of empirical data (see d. and h.). 


a tect 


h. A unique hydrodynamical interpolation technique is introduced which incorporates into 
the theoretical model empirical tidal constants collected from over 2 000 tide-gauge stations 
around the world in a hydrodynamically consistent fashion (see d., f., and g., above). 


i. A new higher order approximation of Arctic Ocean tides is used, that is described in 
Schwiderski (198 1b). 


With these features, the new model was successfully applied to chart the semidiurnal principal 
lunar (M,) ocean tide with the desired accuracy. The technique and accuracy of the model were 
extensively described and discussed in Part | of this report as well as in subsequent publications 
and symposia presentations by the author (Schwiderski 1978a, b; 1979a, b, c, d, e; and 1980). 


The same hydrodynamical interpolation technique has been applied to chart the diurnal 
luni-solar declination ocean (K,) tide with the same relative accuracy as M2. Again, it must be 
emphasized that the enormous accuracy achieved over all open ocean regions diminishes somewhat 
near coastal areas where known empirical data are marginal in quantity and/or quality. 


A complete listing of all sources of empirical ocean tide data, which were interpolated into 
the K, tidal charts, is presented in Appendix A. In the meantime, Section 2 of this report lists 
the significant hydrodynamica! input parameters that specified the constructed K, ocean tide. 
The major features of the global K, tide are discussed in Section 3. A complete numerical display 
is presented in Appendix A where all tidal amplitudes and phases are gridwise tabulated in map- 
like charts. Corange (equi-amplitude) and cotidal (equi-phase) maps of the K, ocean tide are 
plotted in Appendix B. 


tJ 


2. K, OCEAN-TIDE PARAMETERS 


The astronomical diurnal luni-solar declination (K,) equilibrium tide n (or tide-generating 
potential Gn; see Schwiderski, 1978a) at the geographical point (A, ¢) and instant (Y, D, #) is 
determined by 


n = K sin2 @ cos(ot + X + A) (1) 


where 


Il 


9.81 m/sec? earth gravity acceleration 
longitude (east in rad) 

latitude (north in rad) 

(2 1975S) = year number 

day number of year Y(D = | for January 1) 
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universal standard time of day D (in sec) 
= 0.14] 565 m=K, equilibrium tide amplitude 
= 0.72921 - 1074 sec"! =K, tide frequency 
= m(h,, + 90)/180 = K, astronomical argument (in rad) 
= 279,696 68 + 36 000.768 930 4857 + 3.03 + 1074 T? 
© )= mean longitude of the sun relative to Greenwich midnight of day D (in deg) 
T = [27 392.500 528 + 1.000 000 035 6D] /36 525 
D = D+365(Y - 1975) + Int [((Y - 1973)/4] 
Int{x] = integral part of x 
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The corresponding instantaneous ocean partial tide (Schwiderski, 1978a) is determined by 


§ = Ecos(or+X- 6), (2) 


where the local harmonic constants 


wre 


= §&(, 6) = K, ocean tide amplitude (in m) 


and 


fo) 


= &(A. 6) = K, ocean tide Greenwich phase (in rad) 


must be determined, say. by linear interpolation in the tidal charts of Appendix A. 


| 
A simple second-order approximation in the sense of Love and Accad and Pekeris (see Part I. 
| Schwiderski, 1978a, 1979c, and 1980; and Accad and Pekeris, 1978) yields 
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f° = 0.612 and ° = - 0.066078, (3) 


ie., the corresponding terrestrial tide €© and the earth dip §©? (yielding) under the oceanic tidal 
load {, respectively. A more elaborate and probably slightly more accurate earth dip c©? may be 
computed by using Farrell’s Green function (see Farrell, 1972 and 1973; and Schwiderski, 1980). 
In linear superposition, one finds the corresponding instantaneous geocentric partial K, tide: 
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A detailed description of the hydrodynamical-empirical model to compute the ocean tidal 
amplitudes £ and phases 5 (listed in Appendix A) was given in Schwiderski (1978a, 1979c, d, and 
1980). In particular, all model input parameters such as the dimensionless eddy coefficient € (Eq’s. 
103 and 123), the bottom-friction parameter b (Eq’s. 4a and b), and the differencing parameters 
K and K (Eq’s. 64 and 72) were all specified in Schwiderski (1978a) (referenced equations). These 
parameters were determined for M, by extensive trial-and-error computations and remained 
unchanged fer the construction of K,. 


In the computation of the K, tide model, the following mode-dependent parameters were 
used (see referenced equations in Schwiderski, 1978a): 


a. The time step At (Eq’s. 64, 123) 
At = 179.5089 sec (5) 


b. The hydrodynamical interpolation control limits, k,, k2, and k3 (Eq’s. 88, 89. 94, 97, 
and 99) were parily changed to: 


k, = 0.025, k, = 0.040, k; = 0.5. (6) 


It may be noted that the “best” fit to the interpolated empirical tide data required a con- 
siderable reduction of the most significant control parameter k, from the semidiurnal M, and S, 
(see Parts II and III) value of k, = 0.045 to the diurnal K, value of k,; = 0.025. A reduction of this 
sort was anticipated because the empirica! and computed tide data of all semidiurnal tides display 
generally considerably rougher distortions and retardations near boundaries than the diurnal tides. 
Evidently, during the longer time period the diurnal tides are less turbulent and, hence, produce 
less turbulent bottom friction which is reflected in the reduced value of k, (see Part I). 


3. K, OCEAN-TIDE FEATURES 


The entire constructed K, ocean tide is gridwise displayed in map-like amplitude and phase 
tables in Appendix A. The 42° x 71° charts cover the whole globe north of colatitude 169° 
(Antarctica) in three zones: a northern zone N from 0° to 71° colatitude, a middle zone M trom 
48° to 118° colatitude, and a southern zone S from 98° to 168° colatitude. The overlapping 
weographical areas of the tidal charts have been chosen to provide a worldwide coverage for speciul 
applications and to allow the reader to scan the large amplitude and phase charts together in order 
to evaluate their quality and visualize the important tidal features. In addition, a generally super 
ficial overview of some tidal features can be recognized by inspecting the more schematicully 
plotted corange and cotidal maps provided in Appendix B. 


For an easy evaluation of the tidal charts in Appendix A, all hydrodynamically interpolated 
empirical tidal amplitudes and phases have been visibly marked by subbars for all shore data und 
subbrackets for all near-shore deep-sea input constants. Furthermore, the charts display thie 
approximate locations of distant off-shore deep-sea stations by subtides under the computed 
amplitude and phase data. The corresponding empirical data, which were excluded from hydro- 
dynamical interpolation (see Sect. 1 and Schwiderski, 1978a, 1979d, and 1980). are fisted and 
compared with the modeled data in Tables 1, 2, and 3. Finally, the approximate geographical 
locations of the important amphidromic points of zero amplitudes are marked by a circled ®) . 


The tidal charts and maps permit the viewer to follow the tidal waves, that is the high water 
fronts (crests), in forward (or backward) direction, for instance, on their rotation around the 
amphidromic points. In the tidal phase charts of Appendix A, it is best to start from the prom- 
inently visible 0° = 360° or 100° cotidal lines. Since the Greenwich phases specify the time lags 
(in degrees: 15° © 1 hour) of the tidal crests relative to the cresting time of the corresponding 
equilibrium tide along Greenwich meridian, one gathers a vivid impression of the significant global 
and local tidal phenoinena. 


By following the tidal waves on their periodic rotations, one finds these waves passing through 
the specially marked stations in empirically correct time and with the correct height. In fact. all 
over the globe over 2 000 tidal phases and 2 000 amplitudes are coherently integrated. This is 
particularly impressive for the charts of the Pacific Ocean, where the empirical data from so many 
clustered and scattered island stations fit smoothly into the surrounding computed tides. From the 
smoothness features of erratically interpolated tidal data (see Parts I and II), one concludes that 
this result is not an artifact of the interpolation applied but constitutes a vivid manifestation of the 
excellent compatibility of both the empirical and hydrodynamical procedures combined. 


On the basis of this observation, it can again (see Schwiderski, 1978a. b: 1979a, b. d. e: 
1980, and 1981a) be estimated that the K, tidal charts permit a tide prediction with a uniform 
accuracy relative to M, of better than 5 cm anywhere in the open oceans. Naturally, near rough 
ocean basin reliefs (e.g., Arctic and Antarctic shores), where empirical tide (and depth) data are 


marginal in quality and quantity, a somewhat lesser accuracy must be expected. The estimated 
accuracy of the computed K, tide is, of course, fully validated by all 32 empirical tide data from 
distant a shore deep-sea tide gauge staticns, which are listed along with the computed data in 
Tables 1, 2, and 3. The differences (not necessarily errors) range from 0 to 3 cm in amplitudes and 
0° to 15° ( hour) in phases and thus verify the estimated prediction accuracy. 


Table 1. North Atlantic Ocean Deep-Sea Empirical and Modeled K, Tides 


LONGW LATN | EMP MOD ¢ ae EMP5 MODS Aéd | IAPSONR SOURCES 
13°51’ 58°16" / 9 8 -l | 139 137. 2 1.1.37 Cc 
24°43’ =62°50' | 12 10 -~2! 136 138 +2 1.1.29 C 
28°46" 00°12 | 10 «6-1 | 137 141 +4 1.1.30 Cc 
29°58’ 57°01’ | 10 9 -1! 140 142 +2 1.1.31 Cc 
30°10’ 53°39" | 9 7 ~2| 136 139-43 1.1.32 Cc 
25°06" 83°31 6 7 4b} 151 136-15 1.1.33 C 
20°00' ~=—-53°39' ; 8 7 -~1 1 145 132) 413 1.1.34 Cc 
28°11' 48°45’ © 5 =I | 119 122.43 1.1.38 c 
28°09" 45°21" | 5 5 0 | 107 112 +5 1.1.39 C 
27°57" = 44°25" 4 4 0, 94 92 -1; 1.1.40 Cc 
20°05’ 37°09" 5 5 0, 67 65 -2 1.1.41 Cc 
14°15 30°41 / 6 6 O| 61 55-6 1.1.42 Cc 
BAe Wiss ee, Ate i, Seah ce ee es S32 — 
75°38’ 30°47 | 10 10 0} 185 189 4 1.2. 3 C,M 
76°25’ 30°26’ | 10 9 ~1{ 190 195. $5 1.2.11 C,P 
76°48’ 28°27"; 9 9 0} 195 201 +6 1.2.15 Cc 
76°47’ 28°01’ 9 9 0} 196 201. +5 1.2.14 Cc 
67°32" 387 ra 7 <-I | 194 198 +44 1.2. 5 C,Z 
69°45" 28°08’ | 8 7 ~-L | 195 198 43 | «1.2.4 CZ 
69°40’ 27°59’ 8 7 <1} 193 201s +8 1.2. 8 CZ 
69°40' 27°58’ 8 es ae 200," [46.1 29 CZ 
69°20' 26°28" 8 7 Sot | 197 20447 1.2.10 CZ 
69°19' 26°27’ 8 7 -1 | 200 204 #4) 1.2.9 C,Z 


& = Amplitudes (cm) 

6 = Greenwich Phases (deg) 
TAPSO = fnt. Assoc. for the Phys. Soci. of the Oceans 
= Cartwright et al. (£979) 

= Mofjeld (1975) 

= Person (1975) 

= Zetler et al. (1975) 
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Table 2. Northeastern Pacific Ocean Deep-Sea Empirical and Modeled K, Tides 


LONG W LATN |! EMPE MODE AE | EMP56 MOD6 Aé | IAPSO NR SOURCES 
14°29" en 44 44 0 266 268 +2 : 2.1.17 Cc 
135°38' = 53°19"; 43 44 +] 257 259 +2 2.1.16 Cc 
132°47' = 49°35" 43 42 ~! 249 250 +]: 2.4.15 Cc 
145°00' = 34°00" 26 24 —2 224 238 = 86+14 2.1. 9 C1 
145°00’ §=—- 34°00" 27 24 -3 242 238 “4 - I 
124°26' = 27°45’ 30 28 —2 213 213 0 2.1.1 C,M 
129°O}' = 24°47’ 26 25 ~] 222 217 a) 2.1.1 C,M 


= Amplitudes (cm) 
= Greenwich Phases (deg) 
tees = Int. Assoc. for the Phys. Sci. of the Oceans 
= Cartwright et at. (1979) 
= Irish et al. (1971) 
M = Munk et al. (1970) 
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Table 3. Southeast Indian Ocean Deep-Sea Empirical and Modeled K, Tides 


LONG E LATS | EMP MODE AE | EMP& MODS om ess IAPSONR SOURCES 

: i seit ds ice ae 
132°01' = 37°01’ a 17 =l | 231 236 +5 | 4.1. 1 CIS 
132°09' = 50°02’ | - 13 0 232 239 +7 1 4,1. 2 CIS 
132°07' = 60°01’ 18 19 +1 223 231 +8 4.1. 3 CIS 


| 


t= Maieitbities (cm) 

6 = Greenwich Phases (deg) 

[APSO = Int. Assoc. for the Phys. Sci. of the Oceans 
C = Cartwright et al. (1979) 

{S = Irish and Snodgrass (1972) 


Nevertheless, three deep-sea stations with somewhat marginal empirical tide data may be 
pointed out, in order to illustrate the accuracy of the empirical and, hence, computed data. 
For the same Pacific JOSIE II station the original publication by Irish et al. (1971) lists two phase 
values 6 = 224° (Cartwright et al. 1979, [APSO NR 2.1.9) and 6 = 242° for the K, partial tide. 
According to the original paper both values were derived by response analysis from the same 
recorded time sequence. However, the two data were computed with their 18° discrepancy by 
using different references (convolution functions). The inependently computed phase of 6 = 238° 


] 


lies clearly between both empirical values but is distinctly closer to the higher datum. Since the 
lower datum apparently has been preferred by Cartwright et al. (1971), the author attempted to 
enforce the lower phase by hydrodynamical interpolation. The physically unjustified interpolation 
(see Sect. | d. and h.) was abandoned, when the computer experiments generated an almost un- 
etfected surrounding tide and so rejected the lower phase value. 


For the two Atlantic stations with [APSO reference numbers 1.1.33 and 1.1.34 Cartwright 
et al. (1979) lists K, phases which differ by 13° to 15° from the computed phases (see Table 1). 
Though the tides are rather low in amplitude, the mentioned differences are distinctly larger than 
all others of the same area. Therefore, with the same computer experiments mentioned above for 
the Pacific station the author attempted to enforce also the anomalous empirical K, phases by 
hydrodynamical interpolation. These attempts failed also and were abandoned. Now, if one 
follows the almost plane K, tidal wave forward and backward from the distinctly visible 100° 
cotidal (equi-phase) line (see Table 9N in Appendix A), one finds the wave passing properly 
through all other tide gauge stations of this area. Hence, in order to pass also through the two 
stations in question in empirical time, the tide wave must suffer some anomalous retardation in 
the vicinity of these stations. The same anomalous retardation is also directly visible from all 
Atlantic stations on the 53° latitude line listed in Cartwright et al. (1979, IAPSO NR’s 1.1.32 to 
1.1.36). Yet, the depth data at and around the two stations fail to provide a physical cause for the 
empirically determined retardation. Moreover, the tidal waves of the other diurnal and semidiurnal 
constituents except P,, which are all quite similar in shape for this area (see Schwiderski 1979b 
and Cartwright et al. 1979), do not exhibit any symptoms of an analog retardation. The slightly 
anomalous K, and P, phases at the two Atlantic stations can probably be attributed to the distant 
Reykjavik reference station used in the response analysis of these two components with very close 
periods (Schwiderski 1978a). 


From the tidal charts and maps in Appendixes A and B, one concludes that the tidal waves of 
the diurnal K, ocean tide rotate also around amphidromic points qualitatively similar to those 
tidal waves of the semidiurnal species M2 and S, (see Parts If and III). There exist strong and weak 
amphidromes and the pairwise rotations may be compatible or incompatible. However, the dis- 
tribution of the amphidromes displays considerable shifts in position. As was already mentioned 
in Section 2, the distortions and retardations caused by boundary and bottom irregularities are 
generally considerably subdued for the diurnal K, tide when compared to the rougher semidiurna} 
tides as M, and S, in Parts II and III. 


4. CONCLUSIONS 


The hydrodynamical interpolation technique has been applied to construct the diurnal 
luni-solar declination tide (Ky) with a relative accuracy of better than 5 cm anywhere in the open 
oceans. Phe constructed tide is displayed by tabulated charts in Appendix A and by coranye and 
cotidal maps in Appendix B. The major features of the K, tide are discussed in Section 3. A 
comparison with the earlier computed semidiurnal M, and S, tide models reveals qualitative 
similarities. However, the distribution and shape of the amphidromic rotations displays con- 


siderable variations. 
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APPENDIX A 


ATLAS OF 1° x 1° K, OCEAN TIDE AMPLITUDE 
AND PAASE CHARTS FOR 42° x 71° AREAS 


1. GUIDE TO TIDAL CHARTS 


M = m: Longitude Number 

N = n: Colatitude Number 

r = (m- 0.5)°: Geographical Longitude East 

0 = (n- 0.5)°: Geographical Colatitude 

g = &(X,,.0,): Amplitude (in cm) 

ry = 6(\,,. 4): Greenwich Phases (in deg.; 15° = Lh) 


1 i 
on = Amphidromic Points 
= Subbars Mark Empirical Input Data at Shore Stations 
ues = Subbrackets Mark Empirical Input Data at Near-Shore Deep-Sea Stations 
“ = Subtildes Mark Approximately Distant Offshore Deep-Sea Stations with Excluded 


t.mpirical Tide Data Listed in Tables 1, 2, and 3 


tw 


SOURCES OF EMPIRICAL TIDE DATA 
Publications: 


National Ocean Survey (1942), British Admiralty (1977), International Hydrographic Bureau 
(1978), Defant (196)), Miyazaki et al. (1967), Nowroozi et al. (1969), Munk et al. (1970), Zahel 
(1970), Trish et al (1971). Trish and Snodgrass (1972), Nowroozi (1972), Luther and Wunsh 
(1975), Motield (1975), Pearson (1975), Zetler et al. (1975), Cartwright et al. (1979), and Pugh 
(1979). 


Private Communications: 


D.C. Simpson (1977), National Ocean Survey, Rockville, Maryland: S. K. Gill and D. L. Porter 
(197%), National Ocean Survey, Rockville, Maryland: K. Wyrtki (1978), University of Hawaii. 
Honolulu. Hawai, and D. FE. Cartwright and D. Pugh (1978), Institute of Oceanographic Sciences. 
Bidston Observatory, U.K. 
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ATLAS OF CORANGE AND COTIDAL MAPS 
OF THE K, OCEAN TIDE 


Amplitudes & of corange lines in cm. 


Greenwich phases 6 of cotidal lines in 15, 30, 45, 60, 75, 90, 105, 120, 135, 150. 165, 
180, 195, 210, 225, 240, 255, 270, 285, 300, 315, 330, 345, 360 = 0° where 15° = 1| hour. 
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